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High precision orbit determination of Global Positioning System (GPS) satellites is a 
key requirement for GPS-based precise geodetic measurements and precise low-Earth 
orbiter tracking , which are currently being studied at the Jet Propulsion Laboratory 
(JPL). Different strategies for orbit determination have been explored at JPL with data 
from a 1985 GPS field experiment . The most successful strategy uses multi-day arcs for 
orbit determination and includes fine tuning of spacecraft solar pressure coefficients and 
station zenith tropospheric delays using the GPS data. Average rms orbit repeatability 
values for five of the GPS satellites are 1.0, 1.2, and 1.7 m in altitude , cross-track, and 
down-track components when two independent five-day fits are compared. Orbit predic- 
tions up to 24 hours outside the multi-day arcs agree within 4 m of independent solutions 
obtained with well-tracked satellites in the prediction interval. Baseline repeatability 
improves with multi-day as compared to single-day arc orbit solutions. When tropo- 
spheric delay fluctuations are modeled with process noise, significant additional improve- 
ment in baseline repeatability is achieved. For a 246 km baseline, with six-day arc solu- 
tions for GPS orbits, baseline repeatability is 2 parts in 10 8 (0.4-0.6 cm) for east, north, 

and length components and 8 parts in 10 8 for the vertical component. For 1314 and » 

1509 km baselines with the same orbits, baseline repeatability is 2 parts in 10 8 for the 
north components ( 2-3 cm) and 4 parts in 1 0 8 or better for east, length, and vertical 
components. 


\. Introduction 

The NAVSTAR Global Positioning System (GPS), when 
fully operational, will consist of a constellation of 18 satellites 
in 12-hour orbits designed to provide nearly continuous world- 
wide coverage for a variety of civil and military timing and 
positioning applications. For several years, the Jet Propulsion 
Laboratory (JPL) has been actively investigating GPS-based 


measurement systems with the goal of eventually demonstrat- 
ing a sub-meter accuracy capability for GPS orbit determina- 
tion. This capability will provide the basis for centimeter-level 
geodetic studies on a continental scale such as the NASA 
Geodynamics Program’s Caribbean Initiative [1],[2] and 
decimeter-level positioning accuracy for low-Earth orbiters. 
A GPS-based measurement system will ultimately include GPS 
receivers situated at locations throughout the world to provide 
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accurate differential orbit determination for low-Earth orbit- 
ing satellites [3 ] — f 5 ] such as TOPEX, the Earth Observing 
System, the Space Shuttle, and the Space Station, which will 
also be equipped with GPS receivers. Compared to other 
precise positioning systems such as very long baseline inter- 
ferometry (VLBI), satellite laser ranging, and conventional 
ground tracking networks, GPS-based systems offer the poten- 
tial for very high accuracy, an abundance of data, flexibility 
in the placement of ground receivers, and advantageous 
geometry provided by the high altitude and even distribution 
of the GPS satellites. 

This article will discuss strategies for determining precise 
GPS orbits. Different aspects of orbit determination have been 
studied at JPL with 1985 experimental data processed through 
recently developed orbit and baseline estimation software. The 
1985 GPS spring experiment took place between March 29 
and April 5, 1985, representing the first of several GPS field 
tests conducted in 1985 and 1986. More than 20 institutions, 
including JPL, participated in the 1985 spring field test toward 
the goal of testing and evaluating GPS measurement techniques 
and equipment [1]. Data were acquired from seven develop- 
mental GPS satellites in orbit. TI-4100 GPS receivers [6] 
manufactured by Texas Instruments were placed at each of 
the ten ground sites in the continental United States. In addi- 
tion, JPL SERIES-X receivers 1 were operated at the Mojave, 
California, and Owens Valley (OVRO), California, stations, 
and Air Force Geophysical Laboratory (AFGL) receivers, 
forerunners of the Macrometer II [7], were operated at the 
three POLARIS sites (Haystack, Massachusetts, Richmond, 
Florida, and Fort Davis, Texas). The POLARIS sites are 
collocated with VLBI radio telescopes used for geodetic and 
earth orientation studies. Water vapor radiometers (WVRs) 
were available at the Hat Creek, Mojave, and Owens Valley 
sites in California, while at other stations surface meteoro- 
logical information was compiled during the experiment 
and used to correct the data for tropospheric delay. Table 1 
and Fig. 1 summarize the experiment’s configuration and 
the location of the ground receivers. Data were collected for 
about eight hours on each day of the experiment; however, 
because of the geometric constraints of the ground receiver 
configuration, GPS satellites were typically tracked for 2-4 
hours before the receivers switched to new sets of satellites. 

The GPS satellites transmit carrier signals at two L-band 
frequencies (1.22760 and 1.57542 GHz) which are modulated 
by a pseudo-random-noise code (P-code). All receivers in the 


1 R. B. Crow, F. R. Bletzacker, R. J. Najarian, G. H. Purcell, J. I. Statman, 
and J.B. Thomas, “SERIES-X Final Engineering Report,” JPL D-1476 
(internal publication), Jet Propulsion Laboratory, Pasadena, Califor- 
nia, 1984. 


field tests produce “carrier phase” observables from continu- 
ous tracking of the RF carriers. The TI-4100 receivers also 
produce a pseudorange observable by correlating the received 
modulated signal against a local copy of the code. The term 
pseudorange is used since this observable is a measure of the 
light travel time plus clock offsets at the satellite and the 
receiver. The carrier phase observable is analogous to “ambig- 
uous range,” since it measures range biased by an integer 
number of wavelengths. Both AFGL and SERIES-X receivers 
operate without knowledge of the P-code, although the 
SERIES-X receiver also generates a codeless pseudorange 
observable. Data from signals transmitted at the two fre- 
quencies are linearly combined to remove the dominant 
portion of the ionospheric delay, which varies inversely with 
the square of the frequency. Most of the results presented in 
this article are derived from the carrier phase data. Additional 
details on GPS signal structure and characteristics can be 
found in [8] and [9] . 

II. Data Processing 

Data processing at JPL was conducted with the recently 
completed GIPSY (GPS /nferred Positioning SYstem) soft- 
ware. GIPSY includes a comprehensive front end for data 
editing, phase connection, data compression, and atmospheric 
calibrations; PATH-VARY, a program which integrates the 
equations of motion and the variational equations to obtain 
nominal satellite trajectories and transition matrices for 
satellite states and dynamic model parameters; GPSOMC 
(GPS Observed Minus Calculated), a module which calculates 
a very accurate model from best known nominal values and 
computes the pre-fit residuals and measurement partials; a 
U-D factorized batch sequential filter with process noise 
capabilities to perform parameter estimation and covariance 
analyses; and an output processor to display orbits, baselines, 
parameter solutions, and covariances. Clock modeling options 
include explicit single or double differencing, correlated pro- 
cess noise, and quadratic polynomials over specified data arcs. 
Additional details on specific modules in GIPSY can be found 
in [10] . 

Initially, GPS data processing, including parameter estima- 
tion, was performed separately for each day. Then multi-day 
arcs were formed, and longer runs were made covering up to 
six consecutive days. For some days, TI-4100 data were com- 
bined with the AFGL and SERIES-X data for orbit and base- 
line determination. However, since the different receiver types 
were in most cases collocated, there was not much to be gained 
(aside from an averaging effect on data noise) from processing 
all the data together -especially since with multi-day arcs, sys- 
tematic effects as well as data noise become limiting error 
sources. Some of the TI-4100 data were used for determina- 
tion of orbits and baselines independent of those based on the 
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AFGL/SERIES-X network, and these separate solutions were 
compared to assess orbit and baseline repeatability. 


III. Orbit Solution Strategy 

The basic GPS orbit determination strategy includes simul- 
taneous adjustment of the GPS orbits, station and satellite 
clock parameters, selected station locations, zenith tropo- 
spheric delays, range ambiguities (for carrier phase data), and 
solar pressure coefficients. There can be considerable variation 
in the way that these parameters are treated in the filter. In 
this section, various orbit determination strategies are described. 
In the next section, a comparison of orbit and baseline results 
using the different strategies is presented. 

A. The Fiducial Concept 

Fundamental to the GPS orbit determination performed at 
JPL is the fiducial network concept: a network of stations 
whose locations are accurately known from VLBI defines a 
self-consistent coordinate frame to which all GPS orbit and 
baseline solutions are referred. A detailed discussion of geo- 
physical motivations for GPS geodesy, the selection of fidu- 
cial ground sites, and the validation of the fiducial concept 
with GPS measurements from the 1985 spring experiment can 
be found in [10]. Additional discussions of the role of GPS 
orbit determination in high precision geodesy can be found 
in [11] and [12]. An alternative to the fiducial approach, 
the free network approach, is discussed by Beutler et al. [13] . 

The results presented in this article are based on two basic 
fiducial network strategies. The first strategy has four fiducial 
ground stations, typically Haystack, Richmond, Fort Davis, 
and Owens Valley. This provides a well distributed set of 
reference points (Fig. 1). The second strategy has three fidu- 
cial points; the central Fort Davis station is adjusted along 
with the satellite orbits and other estimated parameters. 
Geometry provided by the second strategy is not as effective 
for defining the reference frame. However, in order to test 
baseline repeatability over long (>1000 km) distances, the 
second strategy was used for some of the solution sets. If 
Fort Davis or one of the other fiducial ground site locations 
were susceptible to systematic error or could be improved 
using the GPS data, the second strategy might be expected 
to produce better results. 

B. Solar Radiation Pressure Model 

The earliest GPS orbits determined at JPL showed that 
when fine tuning of the satellite solar pressure coefficients 
was not performed with simultaneous adjustment of station 
locations and GPS states, orbit and baseline repeatability 
was significantly degraded for data arcs longer than eight 
hours (single-day pass). If the solar radiation effects are left 


unmodeled, GPS position errors in down-track can increase to 
over 1 km after 1 -2 weeks of integration of the equations of 
motion. Thus, even a small percentage error in one of the 
solar parameters can result in significant orbit errors. Solar 
radiation pressure is represented in the JPL software with 
the GPS Block I model [14] , sometimes referred to as ROCK4. 
This model implicitly includes spacecraft component shapes, 
orientations, and masses and also accounts for inter-component 
shadowing. In ROCK4, it is assumed that the spacecraft 
remains perfectly oriented with respect to the sun. The space- 
craft-centered coordinates defining the directions for solar 
radiation pressure accelerations (Fig. 2) have the Z-axis posi- 
tive along the antenna directed toward the center of the Earth. 
The Y-axis is along the solar panel support beam, normal to 
the spacecraft-sun direction, and the X-axis is defined relative 
to the other two axes in the sense of a right-handed coordinate 
system. The force model represents the GPS spacecraft with 
13 surfaces, each specified either as a fiat surface or as a 
cylindrical surface. The flat surfaces are defined by length 
and width, while cylindrical surfaces are defined by radius of 
curvature, angle, and length. Each surface’s reflectivity and 
specularity are represented by a number between 0 and 1 to 
account for varying absorption and diffusivity characteristics. 

The solar radiation pressure results in a space vehicle accel- 
eration 


r = P t —— (G a e +Gae) + Ge 

1 2 x x x z z z 7 y y 

- T ps 


(i) 


where e x , e y , e z are unit vectors for the spacecraft-centered 
coordinate system; the scalar P { is a shadow factor, 1 for 
direct sunlight, 0 for umbra, and between 0 and 1 for penum- 
bra; r ps is the spacecraft-sun vector; a x and a z are the space 
vehicle body fixed acceleration components determined by the 
ROCK4 model; k is the nominal Earth-sun distance (1 AU); 
G x and G z are solar pressure coefficient scaling factors; and G y 
is a constant acceleration in the e y direction, in km/sec 2 , often 
referred to as thejy-bias parameter. 


The strategy adopted for fine tuning of the solar pressure 
parameters was to adjust G x , G y , G z per satellite as constant 
parameters over the multi-day arc. Nominal solar radiation 
pressure parameters were taken from the precise post-fit 
ephemeris supplied by the Naval Surface Weapons Center 
(NSWC) [15]. The NSWC ephemeris was considered to be 
accurate in the WGS-72 reference frame to about 15-25 m 
(1 m in altitude and 15 m in each of the cross- and down- 
track components) and was determined from a two-week 
batch fit. The NSWC currently uses a more advanced software 
system known as the Multi-Satellite Filter/Smoother (MSF/S) 
[16] for generating reference GPS ephemerides. The a priori 
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uncertainties for the solar pressure coefficients for fitting at 
JPL from the March 1985 data were assumed to be 25 percent 
in x and z directions and 10“ 12 km/sec 2 for the jp-bias. The 
NSWC fits in March and April 1985 vary from batch to batch 
typically by a few percent in the x and z coefficients and, for 
some satellites, by a significant fraction of 10' 12 km/sec 2 for 
the y- bias. The JPL orbit strategy allowed somewhat more 
freedom in these solar pressure parameters in order to follow 
any shorter term variations which might be overlooked in the 
longer averaging period of the NSWC fits. In addition, the JPL 
orbits were determined with the x and z coefficients estimated 
as independent parameters. Nominally, these two parameters 
have the same value. This extra degree of freedom was left in 
the solution in order to absorb deficiencies in the solar pres- 
sure model and possibly to absorb other long term unmodeled 
accelerations. Some of the known limitations of the ROCK4 
model, discussed by Fliegel et al. [17] , could amount to orbit 
component errors of 4 in or more over a 14-day prediction 
interval. Allowing the JPL solar pressure parameters freedom 
to deviate from the nominal NSWC values was also motivated 
by the need to provide compensation for aliasing due to differ- 
ent parameter estimation strategies at JPL and NSWC. Finally, 
the geopotential model used in the NSWC fits (WGS-72 during 
March 1985) is slightly different from the GEM-L2 used in 
the JPL software. Because of these differences, NSWC and JPL 
force parameter solutions were expected to be slightly different. 

The GPS x and z solar pressure parameter solutions obtained 
from a six-day arc covering March 31 -April 5 were nearly all 
within 10 percent of the NSWC nominal values. The differences 
between the x and z coefficients w'ere also generally less than 
10 percent. The >>-bias adjustments were somewhat larger, 
typically representing 25 percent changes from the nominal 
values. The adjustments for the GPS 6 solar pressure param- 
eters were larger than those for the other satellites by a factor 
of about four. The reasons for this are not known, but these 
relatively large corrections persisted regardless of the manner 
in which the data were analyzed. 

C. Clock Modeling Strategies 

The GIPSY software used at JPL for GPS orbit determina- 
tion offers a number of options for receiver and transmitter 
clock modeling. These options include clocks modeled as poly- 
nomials, explicit single and double differencing, and clock 
behavior modeled as process noise. Most of the results pre- 
sented in this article have all clocks modeled as white noise; at 
each measurement time, the clocks are considered to be inde- 
pendent of their values at other times. The effect is very 
similar to clock elimination through double differencing [18] , 
although the white noise treatment has the advantage of 
eliminating the introduction of correlations between measure- 
ments through explicit differencing. 


D. Treatment of Tropospheric Delays 

The troposphere is represented in GIPSY as a spherical shell 
which adds a delay to an incoming GPS signal: 

p = p z R d (6) + p T RJ6) (2) 

d “w 

where p z is the tropospheric delay at zenith and R is an ana- 
lytic mapping function developed by Lanyi [19] to map 
delays at zenith to the delay at elevation angle 0. The sub- 
scripts d and w refer to the dry and wet components of the 
tropospheric delay, 

A priori values for zenith tropospheric delays were obtained 
from surface meteorology and, where available, from WVR 
measurements. The WVR measurements were converted to 
zenith delays using methods described by Robinson [20] , and 
the surface meteorological measurements were converted to 
zenith delays using the Chao model [21 ] . 

The earliest GPS orbits determined at JPL relied on these 
a priori zenith delays for troposphere compensation. However, 
rms scatter was typically reduced by 50 percent or more when 
wet zenith delay corrections were estimated for each ground 
site in addition to the nominal calibrations. Various approaches 
to the troposphere estimation process were attempted, includ- 
ing not solving for tropospheric parameters at one or more 
WVR sites and varying the a priori uncertainty from 2 to 
20 cm for the wet zenith delay parameter. For multi-day arcs, 
the zenith wet troposphere parameters were initially modeled 
with process noise in such a way that the troposphere values 
from one day to the next were independent but behaved as 
constants on any given day. Subsequently, in an effort to fol- 
low tropospheric fluctuations on time scales ranging from six 
minutes to several hours, the process noise time constraint was 
relaxed so that the tropospheric zenith delay could change 
slowly during the eight hour tracking period. 

Two basic process noise models were used to model tropo- 
spheric fluctuations: colored noise and a random walk. The 
process noise formulation in the GIPSY filter is for first order 
exponentially correlated process noise, as described in detail 
by Bierman [22] and by Wu et al. 2 A brief summary of the 
formulation from those references is given here. 

Let p{t) be the value of a time-varying stochastic parameter 
and let co(f) be a white process noise with zero mean value. 


2 S. C. Wu, W. I. Bertiger, J. S. Border, S. M. Lichten, R. F. Sunseri, 
B. G. Williams, P. J. Wolff, and J. T. Wu, “OASIS Mathematical 
Description,” vol. 10, JPL D-3139 (internal publication), Jet Propul- 
sion Laboratory, Pasadena, California, 1986. 
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The GIPSY filter and smoother are formulated in terms of dis- correlated process noise cases or, for the random walk cases, 

Crete time intervals referred to as batches. At the end of the the cumulative effect of <? dis [Eq. (6)] over one day. 

yth batch, the process noise parameters are updated: 


p. t = m.p. + w. (3) 

*/+ 1 rii v ' 

where Pj ~ p(t.) and p. +1 = p(t. + l ). The exponential multiplier, 
rrij, is defined as 

m. = exp [-(f. + , - f)/r] (4) 

where r is the process noise time constant. 

The discrete time-varying variance for a stochastic param- 
eter is 



IV. Orbit Determination Results 

A. Coordinate Systems 

Nominal orbits used to initialize the filter were supplied by 
the NSWC. In the spring of 1985, these orbits were computed 
in the WGS-72 coordinate system and were expected to be 
accurate to about 15-25 m. The fiducial stations used for the 
precise GPS orbit determination in this experiment are defined 
in the VLBI coordinate system, however, and thus a relatively 
large orbit adjustment was expected which would compensate 
for the difference (mostly a rotation in longitude) between the 
two coordinate systems. An accepted value for the coordinate 
rotation is 0.554 arc sec (2.6 jurad), although published values 
range from 0.5 to 0.8 arc sec [23] . A rotation of 0.554 arc 
sec corresponds (S. C. Wu etal ., see footnote 1) to approxi- 
mately 15 m movement for a ground station and about 70 m 
at GPS altitude. 


= mja 2 ^ + (1 - mj)a 2 ss (5b) 

The steady -state sigma, o ss , is the noise level approached after 
the system has been operating undisturbed for a time much 
greater than r. Together r and o ss define the process noise 
variance in discrete form, q dis . 

The random walk is a special limiting case for Eq. (5). The 
random walk corresponds to r With a random walk, 

there is no steady state and o ss is not bounded, but q d[$ is 
defined in the limiting sense 


<?dis = lim — (6) 


Note that for a random walk, ? dis is related to the Allan 
variance o\ (Ar): 




Table 3 lists various process noise troposphere estimation 
strategies which were attempted in filtering the March 1985 
GPS data. The second column of Table 3 shows the magnitude 
of the process noise for various models— either o ss for the 


In order to empirically determine the coordinate system 
offset with the GPS data, a filter run was set up covering six 
days from March 31 to April 5. The only parameters esti- 
mated in this run were UT1-UTC, X and Y polar motion, 
and white process noise satellite and non-maser station clocks. 
Those station clocks which were using hydrogen masers as 
time standards were modeled in successive runs as linear and 
quadratic polynomials and as uncorrelated white process 
noise. The results were essentially independent of the clock 
model used for the masers. The solution for UT1-UTC was 
about 3.3 jurad, or approximately 0.68 arc sec. The X and Y 
polar motion solution^ were at least an order of magnitude 
smaller. The difference between the JPL GPS result (0.68 arc 
sec) and the published value (0.554 arc sec) can be almost 
totally explained by the difference between the nominal 
earth orientation values supplied by the International Radio 
Interferometric Surveying Subcommission (IRIS) and those 
supplied by the Defense Mapping Agency (DMA). The IRIS 
values are used in GIPSY, while the DMA values were used 
by the NSWC for GPS orbit fitting. The difference between 
the DMA and IRIS UT1-UTC values for March 25, 1985, is 
0.13 arc sec, which, when combined with the 0.554 arc sec 
assumed value for the WGS-72 to VLBI coordinate system 
rotation, is very close to the 0.68 arc sec determined for the 
JPL GPS solution. 

Figure 3 shows two orbit fits for GPS 8 (we use the Navstar 
number to identify the satellites) from the six -day arc. Table 2 
describes the basic orbit determination strategy used for these 
and most subsequent GPS solutions. Nearly all adjusted param- 
eters were estimated with very large a priori uncertainties. The 
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data weights used for all solutions in this article were greater 
than the carrier phase intrinsic receiver noise of several milli- 
meters in order to be consistent with post-fit rms scatter and 
in order to make xl “ 1 , where 



In this expression, n and m are the number of measurements 
and degrees of freedom, z^obs) and z £ .(pred) are the observed 
range and predicted post-fit range for the / th measurement, 
and a- is the measurement noise for the it h measurement. The 
data weights used were based on this criterion and correspond 
to measurement noise of about 0.9-1. 5 cm, depending on the 
length of the data arc. 

The solution plotted in Fig. 3(a) includes the coordinate 
system offset discussed above; most of the rotation shows up 
in the cross- and down-track orbit components. Figure 3(b) 
shows a solution from a filter run identical to the first except 
that a ~3 ^urad coordinate rotation in UTl-UTC (longitude) 
was first included in the model before estimation of param- 
eters. The large cross-track amplitude has been greatly reduced, 
since it is mostly due to the rotation of coordinate systems. 
The down-track run-off can be attributed to a slight overall 
altitude correction to the satellite ephemeris; because the 
orbital period changes with altitude, the down-track correc- 
tion will slowly increase over time. An overall altitude adjust- 
ment of 70 cm, for example, would cause a run-off in the 
down-track component of about 7 m per orbit, roughly what 
is observed in Fig. 3(b). The NSWC nominal orbits are expected 
to be accurate to about 1 meter for the altitude component 
alone, so the results plotted in Fig. 3 are consistent. 

B. Orbit Precision and Accuracy 

Since the effective measurement noise was increased to be 
consistent with the post-fit rms scatter, the formal errors 
from the filter are one measure of the precision of the orbits. 
Figure 4 shows plots of formal orbit errors from a six -day arc 
solution and a single-day arc solution for GPS 8 and GPS 6 fit 
with carrier phase. GPS 8 was well tracked throughout the 
experiment, while GPS 6 was more sparsely tracked, with 
good geometry lasting only a few hours each day. The oscil- 
latory nature of the error magnitude is due to tracking limita- 
tions, since the satellites were observed only from the conti- 
nental United States each day. The single-day run included 
data from nine tracking sites, whereas the six-day run had only 
five to partially compensate for the greater quantity of data in 
the long arc. Figure 4 indicates that orbit precision of about 
0.2-1. 5 m is to be expected from multi-day arcs for well- 
tracked satellites, and about 0.5-3.0m is to be expected for 
more sparsely tracked satellites. For one eight-hour pass 


(single-day arc), the formal errors are somewhat higher- 
1-3 m for GPS 8 and 2-8 m for GPS 6. 

It is important to distinguish between orbit precision and 
accuracy. The formal orbit errors will in general underesti- 
mate orbit accuracy as a result of systematic effects which do 
not appear in post-fit scatter or which cannot be well com- 
pensated for simply by raising the effective measurement 
noise. A number of systematic errors which could grow slowly 
over time might not be apparent from post-fit scatter; such 
systematic effects could result from mismodeling of the earth’s 
geopotential, from long-term spacecraft accelerations due to 
mismodeted solar radiation pressure or other unmodeled 
forces, or from earth orientation errors and fiducial station 
coordinate errors. Thus, comparisons between solutions 
derived from data sets in which different receivers and/or 
different data arcs were used may provide a better measure of 
orbit accuracy than would formal errors. 

In order to better assess GPS orbit precision and accuracy, a 
set of single- and multi-day arcs was set up as shown in Fig. 5. 
Arcs A and B cover three days each and do not overlap in time. 
Arcs C and D cover five days each and do not have any data in 
common, although they overlap in time. Data on March 30 
were also used to obtain GPS orbits completely independently 
of solutions obtained with the six-day arc covering March 31- 
April 5. These different sets of orbits were compared to mea- 
sure orbit repeatability. In all the comparisons involving 
multi-day arcs, the solutions were compared over a “neutral” 
time interval outside the periods during which data were 
taken, so these comparisons provide a rather stringent test of 
the robustness of the orbits and models used to propagate and 
predict up to 24 hours outside the data arcs. 

Figure 6(a) summarizes the GPS orbit repeatability for arcs 
A and B. These results show repeatability based on solutions in 
which the zenith troposphere at each site is adjusted as a con- 
stant bias independent from one day to the next. The overall 
average unweighted rms differences are 1.2, 2.9, and 2.6 m in 
altitude, cross-track, and down-track components for five 
satellites, averaged over a six hour period whose midpoint was 
about eight hours away from the nearest measurements. Fig- 
ure 6(b) shows the corresponding GPS orbit repeatability for 
arcs C and D with constant zenith tropospheric delay param- 
eters estimated. The overall average rms differences are reduced 
to 1.1, 1.5, and 2.4 m in altitude, cross-track, and down-track 
components. Figure 6(c) shows repeatability for arcs C and D 
with tropospheric zenith delay modeled with process noise as 
a random walk (first entry in Table 3). A noticeable improve- 
ment is apparent when Fig. 6(c) is compared to parts (a) and 
(b). With the process noise troposphere strategy, the average 
rms differences are further reduced to 1.0, 1.2, and 1.7 m in 
altitude, cross-track, and down-track components. Figure 7 


6 



shows sample orbit comparisons for GPS 8 and GPS 6 corre- 
sponding to the stochastic troposphere strategy of Fig. 6(c). 

Figures 6 and 7 make a strong case for fine tuning the 
tropospheric delay parameters using process noise. Other 
troposphere process noise models showed similar improve- 
ment over the constant troposphere method, but the random 
walk model for tropospheric delay fluctuations gave the best 
orbit repeatability. Table 3 lists some of the troposphere 
strategies used and notes those which worked best as judged 
by orbit and baseline repeatability. Treuhaft and Lanyi [24] 
have developed a model for tropospheric delay fluctuations 
which behaves in a manner similar to a random walk over time 
scales of several minutes or longer and shows agreement with 
VLBI data. For shorter time scales, the Treuhaft and Lanyi 
model has characteristics of both a Markov process and a 
random walk. However, the tropospheric parameters esti- 
mated from the GPS data are residual delays due to limita- 
tions of the calibrations (from surface meteorology or WVRs) 
whose behavior is not easily modeled. We proceeded on the 
assumption that residual tropospheric delays after calibra- 
tion behave in a manner similar to the tropospheric fluctua- 
tions themselves, although presumably with reduced ampli- 
tudes. All of the rms differences shown in Figs. 6 and 7 are 
consistent with the formal errors (1-a) from the orbit fits. 
Figures 6 and 7 also seem to indicate that with essentially 
similar data and geometry, a five-day arc is preferable to a 
three-day arc for GPS orbit determination. However, it should 
also be noted that arcs A and B have no temporal overlap 
whatsoever; their repeatability is based on orbit prediction 
of up to 12 hours outside the data arcs and should be some- 
what more sensitive to systematic errors which grow with 
time. The C-D repeatability, however, is based on inter- 
polated orbits at times between five and eleven hours from 
measurements but still within the interleaved data arcs. 

Figure 8 shows orbit repeatability using single day arcs. A 
single day typically included 6-7 hours of GPS data; GPS 8 
was tracked for the entire data collection period, while the 
other spacecraft were tracked for periods ranging from one 
hour to four hours. Therefore, in contrast to the multi-day 
arcs, the amount of data and the geometry are very limited in 
a single day arc. The day shown in Fig. 8 is March 30. As 
typically happens during experiments of this kind, there 
were various equipment difficulties on March 30, including 
considerable data loss from the northeast which compromised 
viewing geometry. Although March 30 was not the worst day 
in terms of technical problems, for this experiment its data 
were definitely below average in quality as a consequence of 
data loss. 

Figure 8(a) compares orbits from March 30 data for GPS 8 
in which one orbit solution is from TI receivers and the 
other is from SERIES-X and AFGL receivers. Because one of 


the SERIES-X receivers (Owens Valley) was inoperative on 
March 30, the non-TI data set includes data from only four 
stations. Because of the data loss, the formal errors were rela- 
tively large when station clocks were modeled as white process 
noise. However, when the stations with maser time standards 
had their clocks modeled as polynomials, the formal errors 
decreased to a sufficient degree to make the comparison mean- 
ingful. The AFGL receivers at Richmond and Fort Davis and 
the SERIES-X receiver at Mojave were running on hydrogen 
maser clocks. The Haystack TI receiver used a crystal oscilla- 
tor. The crystal clock fluctuations could be modeled as white 
process noise but not as a linear or quadratic polynomial. 
Figure 8(a) shows agreement at the 1 m level between the TI 
and non-TI orbits for March 30. It is instructive to see that 
even with a compromised data set, the stability of some of 
the station clock standards permitted clock modeling which 
improved the orbit for this satellite. Figure 8(b) compares 
the March 30 GPS 8 orbit with the GPS 8 orbit predicted back 
a full day from the six-day arc covering March 31 -April 5; the 
orbits still agree at the 3 m level. 

Figure 9 shows orbit repeatability with a one-day versus a 
multi-day arc for GPS 6, a satellite representative of those with 
somewhat sparse ground coverage. The single-day/multi-day 
comparison is not as good as with GPS 8, a well-tracked satel- 
lite (see Fig. 8). However, the multi-day/multi-day compari- 
sons for GPS 6 match-and, for some components, even 
slightly surpass-repeatability for GPS 8 (see Figs. 6 and 7). 
Thus, without use of a priori knowledge of satellite trajec- 
tories, 1-2 m orbit repeatability is possible for both well 
tracked and sparsely tracked satellites with multi-day arcs, 
but only well-tracked satellites can achieve this level of repeat- 
ability with a single pass on one day for the limited configura- 
tions of this experiment. This is also reflected in the formal 
errors, which for sparsely tracked GPS were 5-10 m for the 
single day arc but decreased to 1-3 m with the multi-day arcs. 

C. Baseline Repeatability 

Baseline repeatability is often used to assess orbit quality 
in a geodetic context. Baseline repeatability tests have the 
advantage of not requiring any orbit propagation or predic- 
tion. A six -day arc was arranged to test baseline repeatability. 
Five ground receivers (Haystack, Richmond, Fort Davis, 
Mojave, and Owens Valley) were used in these solutions, with 
Fort Davis and Mojave station locations estimated each day 
independently. For each day of the multi-day arc, these sta- 
tion positions were reset with a large a priori uncertainty and 
were estimated anew with the data from that day. One set of 
smoothed orbits was generated based on all the data, including 
data from the two estimated stations. The rms scatters about 
the weighted means of the components of the Mojave-OVRO, 
Mojave-Fort Davis, and Fort Davis-OVRO baselines (246, 
1314, and 1509 km) were calculated. 
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The baseline repeatability results are shown in Fig. 10. Fig- 
ure 10(a) shows the Mojave-OVRO baseline, including a com- 
parison of three different orbit estimation strategies: single- 
day pass solutions, multi-day arc solutions, and multi-day arcs 
with troposphere modeled as process noise. The multi-day arc 
solutions are more consistent than the single-day solutions for 
the 246 km baseline. The stochastic troposphere modeling 
strategy dramatically improves the 246 km baseline repeat- 
ability further; all components are consistent at the 0.5 cm 
level except for the vertical (2.0 cm). Parts (b) and (c), which 
give results for the longer baselines, show less dramatic improve- 
ment in all components than that seen in part (a) for the 
246 km baseline with the troposphere process noise strategy. 
However, the components showing highest rms scatter still 
show moderate improvement. For the shorter baselines, it 
appears that fluctuating tropospheric delay can be a domi- 
nant systematic error source affecting repeatability if left 
unmodeled, but for the continental baselines, orbit errors and 
fiducial network errors may be more important. All baseline 
components over 246-1 509 km show repeatability of 2-4 
parts in 10 8 except for the vertical component for Mojave- 
OVRO (eight parts in 10 8 ). 

It was noticed that Mojave-OVRO repeatability showed a 
slight but significant improvement when the Fort Davis posi- 
tion was adjusted as compared to solutions when Fort Davis 
was fixed as a fiducial station. This should not happen if all 
the fiducial station positions are known a priori to within a 
few centimeters. Therefore, it is possible that both the orbits 
and the baselines being estimated are affected by as yet 
undetermined errors in the Fort Davis or other fiducial refer- 
ence site locations. In general, Fort Davis (either the TI 
receiver or the AFGL receiver) was not adjusted because its 
position was considered well determined from VLB1 tech- 
niques and because it is geometrically advantageous to have 
a fixed reference station centrally located in the United States. 
For the baseline repeatability tests, however. Fort Davis was 
adjusted in order to provide a long (more than 1500 km) 
baseline result. The GPS data should be capable of isolating 
and improving the station(s) with errors in nominal receiver 
coordinates. When this is done, in the future it is likely that 
accuracy for both orbits and baselines will further improve. 

D. Use of the Pseudorange 

Previous covariance analyses [25] have indicated that either 
pseudorange with measurement noise of 10-20 cm (for five 
minute averaging periods) or carrier phase with measurement 
noise of about 1 cm should be capable of providing orbit 
accuracy of 2-3 m with eight hours of data. In the spring 
1985 GPS experiment, the carrier phase noise for six minute 
measurement averages was consistently below 1 cm, and when 
used to determine GPS orbits, those orbits have repeatability 


values slightly better than the 2-3 m predicted from the 
covariance analysis. However, the pseudorange data scatter 
was found to be larger than the expected 20 cm. Figure 11 
shows carrier phase and pseudorange rms post-fit - scatter 
from a solution in which pseudorange and carrier phase 
were processed together. The carrier phase post-fit rms scat- 
ter was uniformly between 0.5 and 0.6 cm for nearly all the 
receivers; however, the pseudorange rms was much larger than 
expected and was highly variable, ranging from 68 cm to as 
high as 225 cm. Ground multipath has been identified as a 
likely source of much of the pseudorange scatter [26] . Satel- 
lite multipath may also have affected the data, but with an 
amplitude of less than 15 cm [27] . 

Tests at JPL showed that a significant fraction of the 
pseudorange scatter could be removed by subtracting a multi- 
path signature from the data; the multipath signature was 
determined by isolating the signal patterns which repeat from 
one day to the next. Nevertheless, the residual scatter was 
still too large for the pseudorange to be useful in a high 
precision orbit and geodetic context. GPS orbits from a 
single-day pass of pseudorange had formal errors (as well as 
repeatability) of about 20-30 m, much higher than the 1-4 m 
(1-2 m for multi-day arcs) levels obtained with carrier phase. 

Despite the high pseudorange data scatter, a limited demon- 
stration was performed to show the potential benefit of pseudo- 
range. With a one-day pass on March 30, the TI-4100 pseudo- 
range and the carrier phase data were processed simultane- 
ously. Each data type was weighted according to its rms scat- 
ter, so nearly all the solution strength came from the high 
quality carrier phase data. Orbit repeatability was measured 
by comparing these combined data type March 30 orbits with 
the carrier-phase-only orbits from a five-station, six-day arc 
covering March 31 -April 5. For satellites with even moderately 
good coverage, only a slight improvement (if any) in orbit 
repeatability was seen with the March 30 combined data type 
orbits compared to March 30 orbits from carrier phase only. 
However, several of the satellites on March 30 had very short 
tracking periods and limited geometric coverage, adverse 
factors which were further exacerbated by data loss which 
occurred for that day. With only carrier phase from a single- 
day pass, these weakly tracked satellites had rather large 
(5-10 m) formal orbit errors. Single-day March 30 orbits for 
these satellites showed significant improvement when pseudo- 
range data was added to the carrier phase measurements, 
despite the high scatter of the pseudorange data. Figure 12 
shows the improvement for GPS 3. Note that the improvement 
is significant for the cross-track and down-track components; 
these components are believed to be weakest with carrier 
phase because of the necessity of solving for range ambiguity 
parameters, which relate the observed range change from car- 
rier phase to the transmitter-receiver absolute range. This 
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weakness in cross-track and down-track is related to a weak- 
ness in baseline solutions for eastern components often seen in 
carrier-phase-only solutions. 

A covariance analysis was performed to evaluate the poten- 
tial of combining high quality pseudorange with carrier phase 
for orbit and baseline determination. The geometry of the 
1985 spring experiment was used for this simulation, with 
seven ground stations, including three fiducial stations. Fig- 
ure 13(a) shows how progressively better pseudorange can 
improve orbits when combined with carrier phase data, partic- 
ularly in the cross-track and down-track components. Adding 
high quality pseudorange to the precise (but ambiguous) car- 
rier phase has an effect similar to that of bias Fixing [18], [28] , 
[29] , since the pseudorange provides an absolute range measure- 
ment which can effectively constrain the carrier phase ambigu- 
ity. Figure 13(b) shows that the east component of a baseline 
improves considerably with the addition of pseudorange to 
carrier phase data if the pseudorange measurement noise is 
less than 50 cm. As a point of comparison, parts (a) and (b) 
also show predicted orbit and baseline errors for carrier range, 
which is carrier phase data for which all range ambiguities 
have been resolved. Note that the mixed data types provide 
accuracies approaching that which would be attainable if 
carrier range were available. With a single eight-hour pass, 
the predicted orbit errors are at the 1 m level or better, and 
the east component of the baseline is determined to about 
3 parts in 10 8 . Since the 1985 spring experiment had a limited 
configuration of ground receivers and satellites, it is expected 
that performance will improve further when the full GPS 
constellation is available. Figure 13 includes an error contri- 
bution from considered parameters based on a 4 cm assumed 
uncertainty in each component of the fiducial station posi- 
tions. The errors from considered parameters are based on 
the sensitivity of the orbits or baselines to unadjusted parame- 
ters; these sensitivities are calculated from the measurement 
partials and geometry [22] . 

If steps are taken in the future to reduce multipath (e.g., 
through antenna design and placement and through use of 
ground plane absorbers), the combined pseudorange and car- 
rier phase data types could become a powerful tool in high 
precision orbit determination and geodetic studies. The accu- 
racies which would be possible with carrier range could be 
approached with combinations of carrier phase and suffi- 
ciently precise pseudorange. 

V. An Orbit Error Budget for GPS 

Formal orbit errors (Fig. 4) are based on the assigned data 
weights and on the a priori covariance for the estimated 
parameters. Since all the estimated parameters (except per- 
haps the tropospheric zenith delays) were basically uncon- 


strained initially in the solutions, the data weights essentially 
determine what the formal orbit errors are. The data weights 
were set to be equal to the measurement noise scaled by a 
factor in accordance with Eq. (8) so that the weights were 
consistent with the post-fit scatter and with the number of 
measurements and estimated parameters. Typically, this 
scaling factor was between 1.0 and 1.5, which indicates that 
some residual systematic errors could be affecting the orbit 
solutions. 

A consider covariance analysis was performed to examine 
the sensitivity of the orbit solutions to unadjusted parameters 
and to various potential systematic errors. As discussed above, 
the errors from considered parameters are combined with the 
formal errors computed in the Filter to produce a total error 
covariance which includes the effects of likely systematic error 
sources whose effects are not compensated for in the basic 
solution strategy. Models for these unestimated parameters 
are used to calculate their effect on the orbits. Systematic 
errors included in this consider analysis were fiducial station 
location errors, errors in Earth orientation parameters, and 
uncertainty in the location of the Earth’s geocenter relative to 
the fiducial stations. The consider analysis covered a six-day 
arc with five ground stations, including fiducials at Owens 
Valley, Richmond, and Haystack. Table 4 describes the 
assumptions for the analysis. The uncertainties in the earth 
orientation parameters specifying UT1-UTC and X and Y 
polar motion are based on the results of Spieth et al. [30] 
and Steppe et al. [31]. These uncertainties apply to earth 
orientation parameters which might be available several 
weeks after a GPS experiment; for real-time applications, 
uncertainties in the earth orientation parameters would be 
expected to be several times higher. Figure 14 shows the 
orbit error breakdown for GPS 8. Note that the expected 
root-sum-square (RSS) orbit errors are close to the l-2m 
orbit repeatabilities’ observed with multi-day arc (Fig. 6[c]). 
The two largest error categories are the computed (formal) 
error and the systematic errors from uncertainties in the 
fiducial station locations. 

The computed error could be reduced by taking more 
measurements, lengthening the arc, including a priori infor- 
mation, or adding pseudorange. The number of measurements 
could be increased by the addition of more ground receivers, 
but as more GPS satellites are launched in the future, the 
amount of data will increase naturally. More a priori infor- 
mation could be included by using satellite and station clock 
models, WVR tropospheric calibrations at more sites, and 
better nominal orbits. Perhaps the most dramatic decrease in 
the computed error could be achieved if bias fixing were possi- 
ble over long baselines [29] . With bias fixing, the east baseline 
component (corresponding to cross- and down-track orbit 
components) computed errors would decrease by at least a 
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factor of two. Even without bias fixing, the addition of pseudo- 
range would similarly reduce the computed error, since it 
tightly constrains the carrier phase ambiguity parameters. 
With some types of receivers, pseudorange is available in addi- 
tion to the carrier phase, but multipath must be reduced con- 
siderably before it is of sufficient quality to obtain sub-meter 
orbits in one pass. Lengthening the data arc beyond six days 
is feasible, but it might be more costly than some of these 
other remedies and might increase sensitivity to systematic 
errors from unmodeled accelerations. 

Even if the computed contribution were reduced by a factor 
of two, the RSS orbit errors would still be over 1 m in cross- 
track and down-track if fiducial station position errors are 
4 cm. Improvement in the fiducial station network requires 
painstaking refinement in the ties between the geodetic monu- 
ments and the GPS and VLBI antennas. Because it appears to 
be one of the limiting error sources, fiducial station accuracy 
is an area of intensive study at JPL. 

The remaining systematic error sources in Fig. 14 are the 
earth orientation parameters (UT1-UTC, X and Y polar 
motion) and the location of the geocenter. These contribute 
relatively little to the total orbit error under the assumptions 
of the consider analysis but could become significant in the 
future as sub-meter GPS accuracy is approached. Note that 
although the earth orientation errors have a small but notice- 
able contribution to the GPS orbit error budget, the effect of 
these errors on differential measurements such as baselines is 
insignificant compared to other error sources [10] . The UT1- 
UTC error contribution in Fig. 14 applies only to the inertial 
reference frame (J2000); in an earth-fixed frame, the UT1- 
UTC error is eliminated in the transformation as long as the 
same value used in the orbit solutions is used to rotate to the 
earth-fixed frame. Most geodetic GPS applications, such as 
baseline determination, are based in an earth-fixed reference 
frame; thus, the UT1-UTC error would have essentially no 
effect, although it would introduce a bias in the GPS orbits. 

Other systematic orbit errors result from unmodeled forces 
and accelerations which can affect the spacecraft trajectories. 
This category of errors includes spacecraft accelerations due to 
gravity mismodeling, gas leaks, solar radiation pressure mis- 


modeling, and atmospheric drag. We believe that unmodeled 
accelerations of this type are small for GPS compared to other 
systematic effects, at least for arcs of one to several days. Pre- 
liminary consider analysis [25] shows very small perturbations 
(less than 5 cm down-track errors) from gravity field mismodel- 
ing for GPS orbits determined from carrier phase. In a future 
study, the analysis of gravity errors will be extended to cover 
longer arcs and to include the gravity field covariance matrix 
to model errors in the gravity coefficients. Since three solar 
radiation pressure coefficients were estimated for each satel- 
lite, we believe that our model has sufficient freedom to com- 
pensate for the solar radiation pressure forces as well as for 
some other minor spacecraft accelerations. When longer data 
arcs are available, these effects will be studied in more detail. 
Finally, the covariance studies found atmospheric drag at GPS 
altitudes to be insignificant, at least for data arcs of up to 
several weeks. 

VI. Summary and Conclusions 

Data from the spring 1985 GPS field test have been pro- 
cessed, and precise GPS orbits have been determined. With 
carrier phase data spanning five days, orbit repeatability is 
l-2m for each component, averaged for five satellites over 
a six hour period during which no data were taken. Baseline 
repeatability over the one-week period for a baseline of 
246 km is about 2 parts in 10 8 (0.4-0.6 cm) for east, north, 
and length components; vertical repeatability is several times 
worse (8 parts in 10 8 ) for this baseline. Repeatability for two 
baselines of 1314 and 1509 km is 2-4 parts in 10 8 (3-5 cm) 
for all components. 

Several refinements to the orbit determination strategies 
were found to be crucial to achieving these levels of repeat- 
ability and accuracy. These include fine tuning the GPS solar 
radiation coefficients and ground station zenith tropospheric 
delays. The time-varying behavior of the troposphere was 
modeled with process noise for the best results. Multi-day arcs 
of three to six days provided better orbits and baselines than 
the eight-hour arcs from single-day passes. A limited demon- 
stration was able to show the potential for further orbit and 
baseline accuracy improvement with combined pseudorange 
and carrier phase data. 
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Table 1. Receiver deployment 


Location 

Receiver 

WVR 

Austin, Texas 

TI4100 

_ 

Dahlgren, Virginia 

TI-4100 

- 

Fort Davis, Texas 

TI4100 

AFGL 

- 

Hat Creek, California 

TI-4100 

Yes 

Haystack, Massachusetts 

TI4100 

AFGL 

- • 

Mammoth Lakes, California 

TI4100 


Mojave, California 

TI4100 
JPL SERIES-X 

Yes 

Owens Valley, California 

TI4100 
JPL SERIES-X 

Yes 

Point Mugu, California 

TI4100 

- 

Richmond, Florida 

TI4100 

AFGL 

- 


Table 2. Parameter estimation strategy 

Parameter 

A priori a 

Satellite positions 

(20, 20, 20) km 

Satellite velocities 

(2, 2, 2) m/s 

Non-fiducial station locations 

(1, 1, 1 ) km 

White noise clocks 

0.1 sec 

Polynomial clocks 

0.1 sec bias 
10' 7 sec/sec rate 
1 0~* 3 sec/sec 2 accel 

Carrier phase bias 

10 sec 

Zenith wet tropospheric delay 

20 cm (no WVR) 
3 cm (WVR) 

Solar radiation pressure 

X, Z, Y-bias coefficients 

25%, 25%, 100% 

Data weights: 1-2 cm 

(white noise clocks all arcs) 

2 cm 

(polynomial clocks 1-day arc) 

Data interval: 

1 carrier 0 meas/360 sec 



Table 3. Stochastic troposphere estimation strategies 



Strategy 

Process noise level 

T 

o p a priori 

1. 

Random walk 

10 cm except Fort Davis/OVRO 
3 cm Fort Davis/OVRO 

oo 

20 cm (SM) 

3 cm (WVR) 

2. 

Random walk 

10 cm except Fort Davis 
3 cm Fort Davis 

oo 

20 cm (SM) 

3 cm (WVR) 

3. 

Random walk 

10 cm except Fort Davis 
3 cm Fort Davis 

oo 

20 cm 

4. 

Random walk 

10 cm 

OO 

20 cm (SM) 

3 cm (WVR) 

5. 

Colored noise 

o ss = 5 cm 

12 hrs (SM) 

24 hrs (WVR) 

20 cm 

6. 

Colored noise 

°ss = 5 cm 

12 hrs 

20 cm 

7. 

Colored noise 

°ss = 10 cm 

96 hrs 

20 cm 

8. 

Random walk 

7.5 cm 

OO 

20 cm 


p: zenith tropospheric delay 
o p : uncertainty in p 

a ss : steady state colored noise sigma Eq. (5) 

Process noise level: 

For random walks: cumulative effect of process noise in Eq. (6)] on o p over 24 hrs 
For colored noise: a ss 

t: colored noise exponential correlation time constant 
WVR: p calibrated with water vapor radiometer measurements 
SM : p calibrated with surface meteorology measurements 

Process noise models listed in order of performance based on baseline/orbit repeatability. 
Models 1-4 performed significantly better than models 5-8 or models with constant p. 


Table 4. Consider analysis assumptions 


Consider parameter 

Consider o 

Fiducial coordinates 

(4, 4, 4) cm 

UT1-UTC 

0.2 msec 

X-pole 

2 masec 

Y-pole 

2 masec 

Geocenter coordinates 

(10,10,10) cm 


Five ground stations, six -day tracking arc 
Fiducials: Owens Valley, Richmond, Haystack 
Station and satellite clock model: white noise 
Troposphere model: random walk 
Carrier 4> data interval: 1 meas/360 sec 
Data weight: 1.5 cm 




■ POLARIS SITES 
• OTHER SITES 


Fig. 1. Locations of GPS receivers during the spring 1985 
GPS experiment 


+x 



-X 


Fig. 2. Local coordinate system used to define GPS solar 
pressure coefficients in the R0CK4 model 


TIME PAST 31 MARCH 1985 04:00:00, hrs 




Fig. 3. Altitude (H), cross-track (C), and down-track (L) orbit 
adjustments for GPS 8 from a six-day arc: (a) with large orbit 
corrections, mostly due to coordinate system offset; (b) with 
~~ 3 jirad coordinate system rotation removed before orbit 
determination 
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POSITION UNCERTAINTY, 


(a) (b) 



TIME PAST 30 MARCH 1985 04:00:00, hrs TIME PAST 30 MARCH 1985 04:00:00, hre 


Fig. 4. Formal orbit uncertainties for satellites tracked from a six-day arc with five ground 
stations and from a single-day (eight-hour) pass with nine ground stations: (a) a well-tracked 
satellite (GPS 8); (b) a sparsely tracked satellite (GPS 6) 
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ALTITUDE CROSS-TRACK DOWN-TRACK 


Fig. 6. Orbit repeatability for five satellites: (a) using arcs A and B 
with constant zenith wet tropospheric delays estimated at each 
station; (b) using arcs C and D with constant zenith wet 
tropospheric delays estimated at each station; and (c) using arcs 
C and D with random walk process noise models for the zenith 
wet tropospheric delays estimated at each station 
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POSITION ADJUSTMENT, 



0 1 2 3 4 5 

TIME PAST2 APRIL 1985 18:00:00, hrs 


Fig. 7. Comparison of orbit solutions for (a) GPS 8 and (b) GPS 6 
mapped to a six-hour period during which no observations were 
taken; these solutions correspond to results plotted in Fig. 6(c) 



0 1 2 3 4 5 


TIME PAST 30 MARCH 1985 04:00:00, hrs 

Fig. 8. Orbit repeatability for GPS 8 with (a) two single-day arcs 
and (b) a single-day arc and multi-day arc solutions 
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TIME PAST 30 MARCH 1985 06:00:00, hrs 


Fig. 9. Comparison of single-day and multi-day orbit solutions 
for GPS 6 



Fig. 10. Baseline repeatability for (a) Mojave-OVRO; 
<b) Mojave-Fort Davis; and (c) OVRO-Fort Davis 




RMS SCATTER, cm 



MOJ FTD HATCR OVR RICH HAY DAHL 



ALTITUDE CROSS-TRACK DOWN-TRACK 


Fig. 12. Improvement in orbit repeatability resulting from 
processing pseudorange data together with the integrated 
doppler from carrier phase for a single-day pass; the single-day 
pass solutions are compared with a six-day, carrier-phase-only 
solution 


Fig. 11. Post-fit rms scatter for (a) carrier phase and (b) 
pseudorange from a solution in which both data types were 
processed together; post-fit scatter is shown for seven TI-4100 
receivers 
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PHASE ONLY +50 cm RANGE +10cm RANGE CARR RANGE 

Fig. 13. Predicted reduction of orbit and baseline errors from 
covariance analysis when carrier phase and pseudorange data 
are processed simultaneously, assuming 4 cm uncertainty in 
fiducial station coordinates: (a) predicted orbit accuracy; 
(b) predicted baseline determination accuracy 



RSS COMPUTED FIDUC UT1-UTC X/Y GEOCENT 

POLE 

Fig. 14. Consider error analysis for carrier phase (GPS 8) showing 
relative contributions to the GPS error budget from data noise 
(computed error), fiducial station coordinate errors, earth 
orientation uncertainty, and geocenter location error; the 
configuration of the spring 1985 GPS experiment was used for 
this analysis 
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